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THE  CRITICAL  EXTERNAL  PRESSURE  OF  CYLINDRICAL 
TUBES  UNDER  UNIFORM  RADIAL  AND  AXIAL  LOAD, 
by  R.  von  Mises, 

Berlin 

Stodola’s  Festschrift,  Zurich,  1929,  pp.  418-430. 

Translated  and  Annotated 
by  D.  F.  Windenborg* 

The  atgnitude  of  buckling  pressure  and  of  the  deformation  of  a  cylindrical 
tube  closed  at  both  ends  and  subjected  to  external  pressure  is  calculated  on  the 
basis  of  elasticity  equations  for  thin  shells.  The  results  of  the  theory  are  com¬ 
pared  with  tests  on  two  series  of  models.  Complete  agreement  is  shown  with  regard 
to  the  lobe  formation;  the  problem  of  determining  the  buckling  pressure  has  not  yet 
been  completely  solved. 


More  then  14  years  ago,  in  view  of  the  conditions  prevailing  with  regard 
to  the  design  cf  fire  tubes  of  boilers,  I  calculated  the  buckling  pressure  of 
cylindrical  tubes  which,  with  fixed  ends,  are  subjected  to  external  pressure.  A 
somewhat  different  problem  is  presented  when  a  circular  cylinder,  closed  at  the 
ends  b?  means  of  bulkheads,  is  deeply  suboerged  in  water  and  thereby  exposed  to 
pressure  from  all  sides,  as  is  approximately  the  case  of  submarines.  The  changes 
in  theory  caused  by  the  introduction  of  axial  load  ire  not  veiy  extensive,  but  the 
altered  numerical  ratios  result  in  the  occurrence  of  phenomena  that  necessarily 
have  a  certain  intrinsic  interest,  especially  concerning  the  large  number  of  waves 
or  bulged  that  becomes  visible  on  the  circumference  of  the  cylinder  at  the  moment 
of  buckling;  see  Fig.  6.  page  12. 

In  the  following,  I  present  a  supplement  to  my  previous  publication^  along 
the  lines  just  indicated.  The  fundamental  equations,  which,  by  the  way,  are  to  be 
found  in  every  text  book  on  the  theory  of  elasticity,  and  the  detailed  calculations 
carried  out  in  the  previous  article,  are  not  repeated;  reference  to  that  article 
is  indicated  cy  the  symbol  Z. 

I  Translator’s  Note:  The  translation  of  the  previous  publication,  of  which  this 
is  the  supplement,  is  given  in  U.  S„  Experimental  Model  Basin  Report  No.  309.  It 
should  be  referred  to  in  c<  ^auction  with  this  article^ 


S.  Emerimontal  Mode]  Basin, 
191*.  p.  750. 


V/unh:lnT .  xi . 


C. 
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1.  DERIVATION  of  the  FUNDAMENTAL  FORMULA. 

Let  us  consider  a  thin-walled,  circular,  cylindrical,  hollo®  body,  closes 
at  the  ends  ty  flat  heads,  and  subjected  to  the  pressure  p  on  its  shell  and  to 

the  compressive  force  P®  on  its 
heads ;  Fig.  1 . 

Let  the  axial  pressure  on  a  unit 
cross-section  of  the  circular  ring  hav¬ 
ing  a  diameter  2a  and  thickness  2h  have 
a  value  p®.  If  the  pressure  on  the 
heads  in  due  to  the  fact  that  the  end 
surfaces  are  subjected  to  the  external 
pressure  p,  it  follows  from: 

FIG.  1  P®  *  IT  a*  p  «  Tr  2a  2h  p® 

that  p®  *  p  a/4h . (1) 

In  order  to  take  into  consideration  the  presence  of  an  end  pressure  p®  the  condi¬ 
tions  of  equilibrium  of  Z  must  be  supplemented.  For  those  equations  are  incom¬ 
plete  inasmuch  as  only  terms  of  the  first  order  in  stress  and  deformation  magni¬ 
tudes  are  retained  in  them.  Products  of  stresses  and  deformations  of  the  second 
order  of  magnitude  are  omitted,  with  the  exception  of  the  products  of  p  by  w  and 
(Z,  Eq.  (9))  since  p,  unlike  all  other  stresses,  is  not  of  the  order  of 
magnitude  of  the  deformations  accompanying  buckling  but  possesses  a  finite  value 
independent  of  them.  We  must  now,  n  an  analogous  manner,  include  also  products 
of  p®  with  deformation  magnitudes  in  our  equations. 


FIG.  2 


Fig.  2  shows  an  element  of  the  shell  in  a  longitudinal  plane  through  the 
axis  of  the  tube,  in  the  deformed  condition.  The  generatrix,  originally  rectilin¬ 
ear  and  parallel  to  the  x-axis,  has  acquired  a  curvature,  which,  except  for  terms 
of  higher  order,  is  measured  by 


Best  Available  Copy 


3 


the  pressure  p*  2h  which  -sots  on  the  surface  of  the  cross-section  of  the  element, 

produces  ft  resultant  in  the  w  direction  (radially  inward,  hence  negative)  equal  to 

-p1  2h  ~j~gr  . . . .  .(2) 

this  expression,  with  opposite  sign,  must  therefore  be  added  on  the  right 
side  of  equation  (9)  in  Z,  'f  we  would  take  into  account  a  general  end  loading. 
Substituting  the  value  of  Eq.  (1)  m  (2)  we  have: 


lie  must,  therefore,  iu 
of  Eq.  (9)  in  Z,  write 

The  same  remarks  appty  to  2q.  (9*)  and  (III)  in  Z.  The  end  loading  has  no  further 
influence  on  the  expressions  for  the  determination  of  the  buckling  pressure,  as 
long  as  we  remain  within  the  limits  of  the  simplifying  assumptions  introduced  in  Z 
or  presupposed. 

It  is  not  difficult  to  follow  the  influence  of  the  supplementary  term, 
equation  (3),  through  to  the  final  formula  in  Z. 

From  Eq.  (11)  of  Z  and  the  abbreviation  a  given  in  Eq.  (12)  of  Z,  we  have 

a*  TF”  *  -«•*"»  <*■  3  . . (4) 

(1  *  effective  length  or  frame  spacing;  see  below,  Section  2). 

Therefore, 

y  (1  -  na  -  «Vz> 

is  to  be  substituted  on  the  right  hand  side  of  Eq.  (Ill*)  of  Z  in  place  of  y Cl  -  n2): 
n  represents  the  number  of  lobes.  This  same  substitution  is  valid  for  the  last 
member  in  the  determinant  (16)  of  Z.  Since  now  y  is  the  only  quantity  in  the  final 
equation  that  contains  p,  an"  therefore,  is  directly  proportional  to  p,  it  follows 
that :  The  buckling  pressure  of  a  tube,  including  the  effect  of  end  load,  may  be 
obtained  by  multiplying  the  pressure  determined  from  Z  by  the  factor 


i  ?*w 

-Pa  a  Yx*~ 


(3) 


tead  of  - 


3aw  . 


in  the  parenthesis  on  the  right  hand  side 


2* w  aa  21 i 

-  ir-T  ip 


na  -  1 
"  pc*  7 
+  it- 1 


.  (55 


statement  applies  primarily  to  the  complete  result  of  the  calculation 
in  Z,  Eq.  (A)  or  (A*),  but  it  may  be  directly  applied  also  to  the 


(8),  since  the  assumptions 


that  led  to  B 


smal Ins 
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x  and  y  defined  by  Eq.  (14)  of  Z - continue  to  hold  generally.  If  we  write  (E) 

somewhat  differently  (by  setting  the  expression  in  the  parenthesis  over  a  denomi¬ 
nator  and  substituting  partially  for  its  value  from  Eq.  (18)  of  Z),  we  get  from 
the  foregoing  expression? 


y  ■  ?T*rr7  M-»") ' 4  a:  ;-|r  —  l> 4  “">•  - 2  * D* 4  * ]■ 


( 


where 


-  i  [l  *  (1  [  2  ■>  (1  y)/oj  =  1  +  1.65/0  +  0.455/0* 

Hi-  ('.  -  rf)jj  +  (i  +  ?.  «-)/>  -  (i  -  Ta)(i  +  frf/Vj 

*  1  +  1.3/w-  1.39  'f*  -  1.417/0*  -f  0.507/O4 


Here  as  in  Z 
h* 


*  *  #■  y *  p  */Jh  •Lir>  (°  ■  ft 


*c» 


and  E  and  tT  are  the  elastic  constants  ((T  »  0.3).  The  determination  of  y  from 
Eq.  (6)  is  not  particularly  f  unbersome  once  the  two  functions  of  p ,  p,  and  p, 
have  been  computed,  and  plac’d  i~.  the  form  of  tables  or  curves  for  fi,  and  fa  . 
Fig.  3  contains  the  two  curves  for  the  entire  range  p  *  0  to  1  sufficiently  ade¬ 
quate  for  all  practical  cases. 

^Translator's  Note;  The  error  noted  in  equation  (B)  of  Z  which  was  carried 
through  to  the  final  formula  (£)  in  Z  has  been  corrected  in  the  derivation  of 
formula  (6).  Hence,  formula  (6)  is  correct  despite  the  error  in  the  original 
equation  (B)  in  zj 
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The  further  amplification  which  was  effected  in  Z  by  the  assumption  of 
only  nail  values  of  p  cannot  be  used  here  since  the  ratio  of  the  radius  to  the 
fraae  spacing  is  sometimes  too  large.  On  the  other  hand,  in  most  calculations  of 
submarine  hulls  we  have  to  do  with  a  larger  number  of  lobes,  n.  It  is  possible, 
therefore,  if  n  equals  at  least  S  or  10,  to  neglect  the  last  two  terms  in  the 
brackets  of  Eq,  (6)  in  comparison  with  the  first  term  (ns  +  ot-a)a.  In  the  same  way, 
we  may  neglect  1  in  the  denominators  in  comparison  with  n8,  by  which  the  error  made 
by  the  previous  assumption  is  partly  compensated.  We  have,  then,  with  sufficient¬ 
ly  close  approximation  for  all  cases  of  practical  importance,  the  final  formula: 


<r*)  ♦  (n*  +*8)8 


(7) 


This  is  even  simpler  than  the  equations  (C)  and  (D)  in  Z.  If  we  substitute  here 
as  before  (T*  0.3,  and  further  for  x  and  y  the  values  from  Eq.  (6*  \  we  obtain 
for  the  required  buckling  pressure  p,  if  2a  designates  the  diameter,  2h  the  wall 
thickness  and  1  the  effective  frame  spacing: 


a 

2  l_ 


2  oc* 

(ns  +  a=T3' 


+  0.73  (n*  +  *8)8  h8/a* 


(8) 


where  <?L  -  TTeJl 

For  n  use  that  7/hole  number  which  makes  the  expression  for  p  a  minimum  for 
a  given  a,  h,  and  1.  E  designates  the  elastic  modulus  of  the  material  as  long 
as  the  stresses  in  the  shell  are  below  the  proportional  limit.  Concerning  the 
applicability  of  the  equations  above  the  proportional  limit,  similar  conditions 
apply  as  in  the  case  of  ordinary  buckling-  of  rods  (aee  below). 

Eq.  (7)  and  (8)  do  not  give  accurhte  values  for  a  tube  of  infinite  length 
ince  here,  as  is  well  known,  the  number  of  lobes,  which  we  have  assumed  to  be 
large,  decreases  to  two. 


2.  DISCUSSION  of  the  FINAL  FORMULA..  TABLES. 

Eq.  (6),  as  well  as  the  simplified  Eq.  (7),  represents  a  straight  line  in 
an  x  -  y  coSrdinate  system^  Each  wave  number,  n,  represents  a  straight  line  for 
a  given  «/ 1,  as  shown  in  Z.  Hence,  y  as  a  function  of  x  will  be  represented  by 
ordinates  of  a  straight  line  polygon.  The  vertexes  of  each  one  of  these  pol¬ 
ygons  increase  in  number  towards  the  origin  of  the  coordinate  system. 

The  significance  of  re  requires  more  detailed  explanation.  According  tc  Z, 

(X,  ~  IT  n/1 . . . . . (?) 

•:  *  r  ■  !;!>;•  r  l<  mshH.j  tat-  her-.-  for  J:  i.nf  t.  ;r,  tr< 
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spacing  of.  two  successive  nodal  points  of  the  lobe  line  appearing  in  a  longitudi¬ 
nal  section  of  the  tubo.  If,  as  is  the  case  in  the  fire  tubes  of  a  boiler,  the 
tube  ends  can  be  regarded  as  fixed  radially,  then  1  can  be  taken  equal  to  the  total 
length  between  the  ends  of  the  tube.  If  the  end  supports  are  yielding,  a  somewhat 
higher  value  of  1  must  be  used.  If  the  tube  is  divided  into  several  belts  by 
sufficiently  rigid  frames,  the  effective  length  1  is  an  average  value  of  the  in¬ 
dividual  frame  spacing.  It  is  possible  to  obtain  more  exact  data  for  determining 

— t 

the  average  value  from  the  known  calculated  results  for  the  buckling  of  rods  with 
several  bays.  However,  since  the  frames  in  submarine  hulls  are  usually  equally 
spaced,  no  difficulty  is  presented. 

It  is  a  question  whether  in  many  cases  a  higher  number  of  waves  between 
two  successive  frames,  that  is  to  say  dividing  1  into  halve3,  thirds,  etc. ,  may 
not  result  in  a  smaller  value  of  the  buckling  pressure.  This  assumption  seems  not 
improbable  when  we  consider  that  the  smallest  possible  number  of  lobes,  n  =  2,  in 
a  circumferential  belt  between  two  frames,  does  not  always  result  in  the  smallest 
critical  pressure.  In  reality  Eq.  (7)  indicates,  in  contrast  to  (C)  and  (D)  of  Z, 
that  for  fixed  n  and  x  the  ordinate  y  does  not  increase  unconditionally  with  cc  . 
For  sufficiently  .small  values  of  x,  y  decreases;  for  example,  if  eC  is  made  three 
times  as  large  while  n  remains  constant.  However,  it  does  not  follow  from  this 
that  the  foregoing  assumption  is  true  since  it  is  not  yet  certain  whether  a  lower 
buckling  pressure  occurs  by  changing  n  while  ct  is  increased.  The  question  to  de¬ 
cide  is  whether  the  polygons  that  are  constructed  for  a  single  OL  intersect  in  the 
region  of  very  small  x  or  whether  the  polygons  always  lie  above  each  other  end 
meet  only  in  the  origin-  The  answer  is  apparent  if  we  consider  that  with  small  x 
on  the  one  hand  the  number  of  lobes  increases  without  limit,  while  on  the  other 
hand  the  vertexes  of  the  polygons  increase  without  limit.  From  the  first  condi¬ 
tion  it  follows  that  Eq.  (7),  by  expanding  in  powers  of  1/n3  and  neglecting  the 
higher  powers,  can  be  simplified  to: 


,  O.t. .CK 


10) 


From  the  second  cordition  it  follows  that  we  may  with  sufficient  accuracy,  in  the 
.  iighborhcod  of  the  origin,  consider  the  polygon  as  the  envelope  of  the  family  of 
straight  lines  determined  by  Eq.  (10)  for  successive  values  of  n.  This  envelope  is 
obtained  through  Eq.  (10)  in  conjunction  with  the  following  equation  obtained  by 
tne  differentiation  of  Eq.  (10)  with  respect  to  n3: 


in 
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Elimination  of  n*  from  Eq,  (10)  and  (11)  gives  the  equation  of  the  envelope: 
y  *  |  («*  +  2)x  *  1.714  ec^x*""  +  |  (<xa  +  2>x  ...  (12) 

As  can  be  seen  from  Eq.  (12),  since  the  ordinates  increase  with  increase  of  &  , 
no  intersection  occurs,  and  we  come  to  the  eonclusion:  The  minisura  buckling 
pressure  always  corresponds  to  a  deformation  without  nodal  points  between  the 
points  of  support  in  the  longitudinal  plane,  and,  therefore,  to  the  smallest  possi¬ 
ble  value  of  oC.  . 

jpTranslators  Note:  Equation  (12)  derived  to  justify  the  conclusion  in  the  preced¬ 
ing  paragraph  is  a  good  approximation  to  the  envelope  of  Eq.  (7)  in  the  neighbor¬ 
hood  of  the  origin,  that  is,  for  small  values  of  x  and  ec.  The  exact  equation  of 
the  envelope  ha3  a  very  important  property:  It  can  be  used  in  place  of  the  original 
equation  (7)  to  determine  y  with  sufficient  accuracy  for  any  given  values  of  x 
and  X-.  Graphically,  this  fact  is  evident  from  the  proximity  of  the  envelope  and 
those  parts  of  the  family  of  straight  lines  used  to  determine  y  represented  by  Eq. 

(7)  with  n  as  a  parameter.  Analytically  it  is  evident  from  the  fact  that,  the 
equation  of  the  envelope,  obtained  by  eliminating  n  between  Eq.  (7)  and  (11),  is 
in  effect  Eq.  (7)  with  that  value  of  n  substituted  which  for  any  given  x  and 
makes  y  a  minimum. 

Eq.  (12)  is  an  approximate  equation  of  the  envelope  of  Eq.  (7),  valid  only 
in  the  neighborhood  of  the  origin,  because  of  the  nature  of  the  assumptions  which 
led  to  Eq.  (10).  The  equation  of  the  envelope  of  Eq.  (7)  should  be  derived  directly 
from  Eq,  (7)  itself.  When  so  derived,  its  usefulness  is  no  longer  limited  to  the 
neighborhood  of  the  origin.  This  equation  can  be  very  accurately  derived  as  follows: 

Differentiating  Eq.  (7)  with  respect  to  either  n,  (n*  +  <£3),  or  p  and 
equating  to  zero  we  get 

(n*  +  X  -  (n*  +  ccsf  <X2x  -  3  (n2  +  <*3)  <*4<1  -  rs)  +  -  <7a)  =  0  (;1a) 

* 

The  solution  cf  (11a)  for  n  gives  that  value  of  n  which  will  make  Eq.  (7)  a 
minimum  for  any  given  oc  and  x.  A  very  approximate  solution  can  be  readily  ob¬ 
tained.  Rewriting  Eq.  (11a) 

(n:i  +  £(n2  +  <xs)  x  -  <£3x]  ~  3^*0  -  G""3)  f  n3  +  ct 2  _  <Xs/3j  =  0 

(n>  +  <*S}4  »  3  -  r2](n*_±  2XrH/j) 

nx 
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la  the  majority  of  practical  cases,  the  value  of  the  ratio  #/n  lies  between 
2/3  and  1/3.  Placing  >*  1/2  in  Eq.  (11c)  and  substituting  (11c)  in  (12a)  we  ob¬ 


tain  the  siaple  expression 

,  .  J,7M_gJfc  . (a) 

1  -  0.374  «tfT 

lith  «/n  •  2/3  v«  obtain  the  slightly  altered  equation 

.  1.734  &JE1 . (b) 

1  -  0.365  «tvx 

* 

while  <x/ n  *  1/3  gives 

.  i-7i?  *  y? . (c, 

1  -  0.332 a: yx 


Any  number  of  similar  equations  may  be  obtained  ty  using  different  values  of 
the  ratio  <*/n.  As  a  limiting  case  ocj n  *  0  gives  the  equation 


1.714  ec  v'x*' 
1  -  0.389*f'T 


•  •  «  .  td  ) 


which  differs  but  slightly  from  (c). 

If  equation  (d)  be  expanded  in  a  series,  the  first  two  terms  are  identical 
with  the  first  two  terms  of  equation  (12).  Moreover,  the  expansion  of  any  of  the 
equations  (a)  to  (d)  yields  a  series  all  terms  of  which  are  positive.  Therefore, 
the  same  conclusion  can  be  drawn  from  these  equations  as  from  Eq.  (12)  concerning 
the  impossibility  of  nodal  points  between  points  of  support  in  the  longitudinal 
plane.  It  ia  evident  tint.  Eq,  (a)  -  (c)  are  such  better  approximations  to  the 
envelope  than  Eq.  (12).  Each  gives  extremely  accurate  values  in  that  particular 
region  of  «c/n  for  which  it  was  developed  and  good  values  for  a  considerable  rang© 
of  i^/n. 

To  determine  just  how  closely  these  approximate  equations  oheok  with  Eq. 

(6),  about  70  values  of  y,  covering  the  region  x  »  4/3.10"4  to  12. IQ"4  (t/d  *  0,002 
to  0.006)  and  <*.  *  1/4  to  16  (1/d  *  6.28  to  0.1),  were  computed  by  each  of  the*  sis 
>-v=  it.1 (/).  -C7)p  (a),  (b),  (c),  and  (d).  These  values  and  their  percentage 

• «,  -  %  •  •'  viiluHS  computed  by  So.  (6)  are  given  in  Table®  II  and  HI, 
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Of  the  four  approximate  equations  tried,  equation  (a)  shows  the  best  agree¬ 
ment  with  the  exact  equation  (6).  The  agreement  is  excellent  for  all  values  of 
1/4  less  than  2.  The  mean  of  the  absolute  deviations  throughout  the  range  of  1/d  = 
0.1  to  2  is  0.5  per  cent.  The  maximum  discrepancy  is  less  than  3  per  cent  for 
1/d  *  2  and  is  considerably  below  2  per  cent  for  all  values  of  1/d  equal  to  or  less 
than  1.6  and  is  less  than  1  per  cent  for  all  values  of  1/d  less  than  1.  In  fact, 
Eq.  (a)  gives  better  agreement  with  Eq.  (6)  than  does  Eq.  (7)  from  which  it  was 
derived! 

Equation  (a)  can  be  put  in  a  more  convenient  form  by  substituting  values  of 
y,  x,  and  oc  from  Eq.  (61)  and  (9),  whence 


2.424  E 

:i>  <?*) 


(h/a  ? 


or  (0*=  0.3) 


-  0.447  Vh7a 
,  2.60  E  v  , 


.  .(a') 


r  j  -  0.45  Vt/d 

Again  it  might  be  pointed  out  that  the  usefulness  of  Eq.  (a1)  is  not  con¬ 
fined  to  the  neighborhood  of  the  origin.  It  is  a  simple  equation,  independent  of 
n,  which  determines  p  with  a  high  degree  of  accuracy  for  a  large  range  of  values 
of  1/d  and  t/d.  Equation  (a1)  may  replace  equation  (6)  in  all  practical  computations 
where  1/d  is  not  greater  than  2.  When  1/d  is  less  than  1  the  discrepancy  will  be 
less  than  1  per  centQ 

In  order  to  facilitate  the  practical  application  of  the  results  derived 
above,  one  may  use  to  advantage  a  graph,  Fig.  4,  that,  contains  the  polygons  (they 
appear  as  curves)  for  the  values  OC  *  2,  4,  6,  .  .  .  ...  to  20,  within  the 
range  u?  to  x  3  6  x  10'4  corresponding  to  the  ratio  a/1  from  approximately  2/3  to 
6  (1/d  3  0.3  to  0.08)  and  to  the  ratio  h/a  up  to  about  0.0045.  The  straight  lines 
are  computed  throughout  by  the  complete  equation  (6).  In  the  neighborhood  of  the 
origin,  the  polygons  are  made  up  of  curves  determined  by  Eq.  (12).  The  corre¬ 
sponding  angles  of  the  polygons  are  connected  by  dotted  lines,  forming  quadrangular 
areas  for  the  various  wave  numbers.  The  separate  vertical  straight  lines  correspond 
to  the  constant  values  c.f  1900  h/a  written  beneath;  the  values  of  the  ratios 
ifa/l  are  written  on  the  3ides  of  the  polygons.  On  the  axis  of  ordinates,  a  second 
scale  is  represented  beside  the  y-snale  which,  under  the  assumption  ox’  an  el art- r 
fivdu  !*ir.  K  2  2, 125, 000  kg/cma  (30.22  x  10*  lb.  per  sq.  in.),  gives  the  value  of 


--  2.34  x  10  y 


i/l,  '  ’-2  !  ,'  T  so.  j.i.  )  The  values  in  t.'n 
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buckling  of  rods.  According  io  these,  as  is  trail  known,  the  critical  pressure  for 
a  rod  of  slenderness  ratio  r/1  is  determined  by  means  of  two  constants  e  and  K  by 
the  formula 

k  *  K(1  -  c  1/r)  . . . . . . (14) 

where  for  mild  steel  K  =  3100  kg/ cm®  (44,000  lb.  oer  sq.  in.)  and  c  s  0.00368. 
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FIG.  4.  GRAPH  FOR  DETERMINATION  OF  BUCKLING  PRESSURE. 
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FIG.  5.  Relationship 
between  K  and  y. 
which  the  buckling  pressure  p  is  obtained. 


were  to  be  substituted  in  the  otherwise  unaltered 
Euler  buckling  formula.  We  have  employed  this 
expression  for  E  in  Eq.  (13)  for  the  conversion 
of  y  into  k  for  k  ^  1900,  y  -  0.090309,  and 
thereby  obtained  the  k-scale  for  the  unelastic 
region.  The  relation  between  y  and  k  is  repre¬ 
sented  in  Fig.  5* 

The  application  of  the  graphs  ia  as  follows. 
For  the  given  value  of  the  ratio,  wall  thickness 
to  diameter,  i.e.  h/a,  we  seek  the  corresponding 
ordinate,  as  well  as  the  corresponding  polygon, 
for  the  given  ratio  of  the  radius  to  effective 
length  a/1.  The  ordinate  of  the  intersection  of 
these  two  lines  gives  the  value  of  y  or  of  k  from 


or  p 


2h 

a 

■ 


k 


(16) 


If  no  lines  exist  for  the  exact  ratio  given,  y  can  be  determined  by  interpolation. 
In  any  event,  the  graph  indicates  the  number  of  lobes,  n,  which  is  determinHtive 
for  the  given  case.  For  a  known  n,  it  is  then  not  difficult  to  compute  the  exact 
value  of  y  by  equations  (6),  (7),  or  (8). 

The  following  numerical  table  was  obtained  by  readings  from  the  graph.  It 
is  ’ alid  both  below  and  above  the  proportional  limit,  providing  the  cylindrical 
shell  is  manufactured  from  medium  steel  plate. 


Table  I. 


Critical  pressure  p  in  kg/cm*  for  various  wall-thickness  ratios, 
h/a,  and  length  ratios,  s/1. 


a/1 

1000  ’n/a 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

2 

“ 

0.637 

0.045 

0.23 

0.6 

1.3 

2.2 

3.6 

5.3 

7.4 

4 

1.273 

0.080 

0.46 

1.3 

2.6 

4.6 

7.4 

10.9 

15.3 

6 

1.010 

0.125 

0.7? 

2.0 

4.0 

7.2 

Vi.  5 

14.2 

17.2 

# 

2.  id/ 

0,175 

1.0 

2, 7 

’.6 

9.6 

12.5 

15.4 

18.3 

10 

?.  I  S? 

0.22 

1.2 

3.5 

7.2 

10.3 

13.2 

16.1 

1  ;■ 

3.800 

0. 26 

i.; 

4.2 

8.0 

10.9 

15.8 

\A 

/. ,  .1  f  f  -. 

o. 

1 , 8 

5.1 

8.4 

11.2 

V(-| 

/.  < 

_  pf 

8.v 

11.6 
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If  we  proceed  purely  analytically,  without  the  use  of  tables,  we  must  use 
the  equation 

. 

in  the  unelastic  region,  which  follows  from  the  substitution  of  Eq.  (15)  in 
Eq.  (13)  and  the  elimination  of  k.  The  velue  of  y  must  here  be  deteimired  by  the 
equations  in  Section  1.  The  numerical  values  in  Eq.  (17)  are  for  medium  steel. 


(17) 


3.  COMPARISON  with  EXPERIMENTAL  RESULTS. 

Two  series  of  experiments  were  carried  out  in  1918  with  large  medium  steel 
tubes  welded  along  the  longitudinal  seam,  and  subjected  to  external  water  pressure. 
The  first  series  included  tubes  with  the , wall  thickness-ratio 

h/a  «  1/400  »  0.0025;  x  »  hV3«*  3  2.083  »  10“ 4 

and  the  three  length  ratios 

a/1  =  400/120,  400/180,  400/240  whence, 

*  =ra/l  3  10.47,  6.98,  5.24 


le  obtain  from  the  graph,  Fig.  4,  at  the  ordinate  h/a  *  0.0025,  the  values  of 
the  number  of  lobes,  n,  in  the  three  cases: 

(a)  n  3  16  (b)  n  »  14  (c)  n  *  13 

[Translator's  Note:  For  case  (c),  n  lies  so  close  to  the  border  line  between  12 
and  13  that  practically  it  might  be  either.  (Actually,  n  3  12  is  determinative?] 
The  experimental  number  of  lobes  from  the  tests  on  four  different  sizes  of 
models  of  diametere  800-,  1200,  1600,  2400  mm  (31.5  in.,  47.2  in.,  63.0  in., 

94.5  in.)  are 


Modal  Size 

Case  (a) 

(b) 

(c) 

I 

(D  3  800  mm) 

n  =  17 

n  3  15 

n  3  14 

II 

(D  =  1200  ran) 

17 

15 

13  -  14 

III 

(D  =  1600  mm) 

— 

14  - 

15 

— 

IV 

(D  3  2400  mm) 

16 

14  - 

15 

— 

tion 

model 

■■f 


it 


We  may  conclude,  therefore,  that  the  agreement  is  almost  exact.  In  addi- 
is  worthy  of  note  that  the  agreement  appears  to  improve  for  the  larger 
The  differences  between  the  number  of  .o’^es  observed  in  the  various  sizes 
Is,  » a  well  os  the  departures  from  the  calculated  number,  hardly  exceed  the 


•f  naav.>:  *  ■  >/ .■‘H mental  error. 


Fig.  6  shows  the  buckled  model  in 
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PIG.  6.  A  Portion  of  the  Experimental  Tube  IVa  with  1o  lobes 

Around  the  Circumference. 

The  situation  is  quite  different  in  determining  the  buckling  pressure  it¬ 
self.  Here  the  experimental  values  obtained  for  the  individual  sizes  of  models 
differ  very  materially  from  each  other.  Since  the  collapsing  pressures  increase 
directly  as  the  size  of  the  model,  it  might  be  expected  that  the  calculated  values 
should  correspond  to  a  limiting  case  of  an  infinitely  large  test  body.  According 
to  Eq.  (6),  using  the  approximation  of  Jig.  3,  the  computations  show: 

[Translator’s  Note:  Numerous  minor  numerical  errors,  probably  due  to  the  use 
of  a  slide  rule,  have  been  corrected  by  the  translator  without  specific  reference 
in  each  case.  Changes  from  the  original  text  are  indicated  *  J 


Case  (a)  oC  =  10.47*  n  »  16 


f>  = 


gca  _ 
H'*  +  ns 


0.30 


ft  =  1.53 


0.91  x  0.09 


ft=  1.23 

+  '365.6s*  -  2  x  1.53  x  256  +  1.23)  x  10"6 


y  =  0.000264*  C.  000894*  =  0.00116 

: ien  above  the  limit  0. 00081 ,  p  must  be  determined  from  v  by 
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^Translator' s  Tote:  Substituting  the  above  values  CC  *  10.47,  x  -  2.083  x  1(f* 
in  (a)  gives 

.  1.718  x  10.47  «  (2.033  x  IQ-*)*  .  0.OOII6 

1  -  0.374  x  10.47  x  (2.053  x  10~*r 

Case  (b)  cc-  6.98*  '  n  *  14  p  =  0.199  3  1.35  1.20 

y  »  M3^222^  (244.7*  -  2  x  1.35  x  196  +  1.20)  x  10~6 

y  =  0.000164  ♦  0.000564  *  0.00073 

This  value  lies  below  the  limit- and  therefore  according  to  Eq.  (13) 

p  =*  0.005  x  0.00073  -  0.5  kg/cm*  »  8.2*  atm. 

[Translator's  Note:  Substituting  E  -  2.125  x  10*  kg.  per  sq.  cm. 
t/d  =  0.0025  1/d  =  2'  x^00  *  0'?25  directly  in  Eq.  (a«) 


2.60„x, 2.125  xJ0<_x^0.0Q25j^3  3.53  kg/c*  ,  g.26  at^j 
0.225  -  0.45  x  y  0.0025  -J 


Case  (c). 


5.24,  n  *  13,  0.140,  /*,  =  1.24*,  ^=1.2 


y  »  P.^1.,,^^0196  +  2^  (196.5**  -  2  x  1.24*  x  169  +  1.2)  x  10** 
»  0.000098  +  0.000438  *  0.000536. 

0.005 x  0. OOj;36  *  6.26*  kg,  per  sq.cm.  »  6.06*  atn. 


(Translator's  Note:  The  use  of  n  =  12  instead  of  n  *  13  in  the  above  problem 
gives  a  slightly  smaller  value  of  y. 

t*  =  5.24,  r.  =  12,  f>  =  0.160,  f,  =  1.28  flt=  1.2 
y  "  q'-^V^T2^  +  (171‘5*  “  2  x  x  144  +  1.2)  x  10~6 


*  0.000149  ♦  0.000386  =  0.000535 

However ,  ir.  practice,  collapse  may  occur  in  either  12  or  13  lobes. 

-  { i;  7  )  give? 


b.  22  kg.  r>‘?-r 


-  6.02  a 
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Theoretical  Values:  (a)  p  3  10.2*  (b)  p  =  8.2*  .c)  d  -  . 

Experimental  Values:  (a)  (b)  ^ 

Model  Sire  I  6.1  3.9  3-35 

II  7.0  5.1  4*15 

III  -  6*°  - 

IV  9.5  7- 

(The  entirely  disproportionate  value  of  Ilia 

is  omitted  here.  Fig.  7  gives  a  general  idea 

of  the  distribution  of  these  values.  Here  we 

see  that  the  calculated  values  may  very  veil 

_ be  regarded  as  obtainable  asymptotic  limits. 

yftr"  There  remains  only  the  question  as  to  what 

S'  -■ ~~ '  causes  the  smaller  models  to  collapse  so  far 

’^-^z*3*** -  l3elow  the  theoretlcal  collapsing  pressure.* 

It  seems  most  probable  that  the  welding  of 

the  thin  plating  disturbs  the  symmetry  of  the 

circular  form,  or  the  homogeneity  of  the  ma- 

**  terial  itself.  It  is  certain  that  in  the 

___ __ _ _ .  - _ _ _  case  of  buckling  with  higher  number  of  lobes, 

FIG.  7.  Summary  of  Experimental  even  very  slight  irregularities  may  have  a 
Data. 

decisive  effect. 

The  second  group  of  experiments  concerns  tubes  with  the  wall  thickness-ratio 
h/a  *  0.00406,  x  a  h*/3a*  s  5.49**  10  6 

and  wit-j  spacing  of  400,  500,  and  600  mm  (15.7  in.,  19.7  in.,  23.6  ir.. )  and 
a  tube  radius  of  800  mm  (31.5  in.).  In  both  the  first  two  cases  the  frames 
failed  under  the  test.  In  the  third  caB£,  likewise,  it  must  be  assumed  as  certain, 
because  of  the  magnitude  of  the  observed  buckling  pressure,  that  the  limit  of 
strength  was  reached  through  failure  of  the  frame. 

[Translator's  Note:  It  is  understood  that  the  testing  arrangements  were  such 
as  to  make  it  impossible  to  observe  the  inside  of  the  model  during  the  testsTj 

As  w*.  know  from  a  theory  which  will  not  be  further  discussed  here,  if  the 
frames  arc  too  weak,  the  buckling  pressure  must  correspond  to  that  of  a  tube  vith 
an  effective-  length  of  double  the  frame  spacing  or  even  higher.  We  can,  therefore, 
test  the  reliability  of  our  theory  by  computing  the  buckling  pressure  for  tubes  of 
the  length  1  *  800,  1000  and  1200  min  and  -cupare  them  with  the  test  results,  with 
help  of  the  graph  we  find: 

(1)  <5t  =  '3.i4  (2)  <£=2.51  (3)  CL  =>2.09 

n  8  n  -  8  n  =  7 


FIG.  7.  Summary  of  Experimental 
Data. 
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y  „  0121^_0|018.+  (73.9*  -  2  x  1.24  x  64  +  1.2)  x  10  6 

»  0.000241  +  0.000428*  =  0.000669* 

p  .  0.00812  x  *  0.000869*  -  12.7*  **■  «!•“•  =  12'3*  ata 

/f«wl.tor'«  Note:  This  valu.  of  p  is  too  hieh  since  n  «  9  is  determinative 
instead  of  n  *  8.  0sin«  n  •  9,  /•*  0.109.  /1i='-19.  /*»"  1-1 

y  .  O.^xO.OHS  ,  (90.9*  -  2  x  1.19  x  81  ♦  1.1)  x  10T* 

0.000127  *  0.000522  =  0.000649 

p  -  0.00812  X  2,1  *  O.OOOM9  -  12.31  kg.  per  sq.cm.  »  11.9  »*■»• 


The  use  of  equation  (a')  gives 

v  o  IS*  t  Iff  r.  (0.00406)**  a  12  ?1  per  sq.  Co.  -  11.9  atm. 

P  "  0.5  -  0.45  y  0.00406 

It  will  be  noted  that  Eq.  (a*)  gives  the  correct  collapsing  pressure  independent 
,f  the  number  of  lobes,  while  the  use  of  Eq.  (6)  with  an  incorrect  value  of  n 
g^vec  a  collapsing  pressure  which  is  too  high^ 

(?>  f  »  0.0896*  ^,=  1.16  f*i=1'1 

y  Q-91^^0080^  *  (70.3**  -  2  x  l.lo  x  64  +  1.D  x  IQ’6 


=  0.000110*  +  0.000398*  *  0.000508* 

p  -  0.00812  x  x  S-000503*  =  pei‘  6q-  ^  =  9‘3*  Stm 

[Translator's  Note:  The  use  of  Eq.  (af)  gives: 

_  2.60  x  2.125  x  10fc  x  /g*$4^-£  =»  9,73  kg.  per  sq.cm.  *  9.4  atm. 

P - 0.625  -  614?  y0.0°40b  6  -J 

(3)  0.0820  1-13  /*»,-  1*1 

y  ,  0,^1^^0067  +  1^  (53.4*  -  2  x  1  13  x  49  +  1.1)  x  10"6 

M  ryrnl/l  *'  0.000300*  =  0.000421* 
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[Translator1 a  Note:  The  use  of  Eq.  (a’)  gives: 

p  a  *.-2.0 £lJL12. — =  8.04  kg.  per  sq.cis.  =  7.?  atm?) 

0.750  -  0.45  x  y  0.00406  J 

The  experimental  pressures  12.2,  9.2,  and  8.5  atmospheres  are  to  be  compared 
with  the  three  calculated  values  thus  obtained,  11.9,  9.?,  and  7.7  atmospheres. 

[Translator’s  Note:  These  values  are  given  as  12.8,  9.5,  and  8.1  in  the  textTJ 
fle  see  that  there  is  good  agreement. 

The  n»imber  of  lobes,  n  =  14  —  15,  observed  in  the  third  test,  as  compared  to 
the  theoretical  number,  n  =  7,  cannot  be  explained.  Probably  it  is  due  to  an  error 
made  by  the  observer.  In  both  the  other  tests  it  was  impossible  to  determine  the 
number  of  lobes  after  failure. 

Conclusions  (added  fcy  translator). 

Equation  (6)  gives  a  solution  for  the  collapsing  pressure  of  a  circular 
cylindrical  vessel  closed  at  the  ends  by  flat  heads  and  subjected  to  uniform  hy¬ 
drostatic  pressure  on  both  shell  and  heads.  As  will  be  observed,  the  determination 
of  the  correct  collapsing  pressure  depends  not  only  upon  the  length,  diameter  and 
thickness,  but  also  upon  the  number  of  lobes,  n,  into  which  a  circumferential  belt 
of  the  vessel,  between  stiffening  rings,  divides  itself  at  collapse.  Equation  (6) 
gives  a  different  value  of  p  for  each  assumed  value  of  n  and  collapse  occurs  in 
that  number  of  lobes  for  which  p  is  a  minimum.  The  use  of  Eq.  (6),  therefore, 
involves  the  determination  of  the  value  of  n  which  gives  a  minimum  value  of  p, 
either  by  actual  substitution  in  that  equation  or  by  the  use  of  curves  previously 
prepared,  such  as  given  in  Fig.  4.  However,  the  determination  of  n  by  Fig.  4  is 
not  always  accurate,  as  was  shown  in  the  solution  of  example  1,  page  15,  where  the 
author  by  using  n  =  8  instead  of  the  correct  value  n  =  9  obtained  a  collapsing 
pressure  which  was  considerably  too  high.  Fig.  8  on  page  18  is  drawn  to  a  differ¬ 
ent  scale  and  gives  a  very  accurate  determination  of  n. 

Equation  (7)  is  simpler  than  Eq.  (6)  and  the  collapsing  pressures  calculated 
by  it  differ  from  those  calculated  by  Eq.  (6)  by  less  than  1  per  cent  for  all  values 
of  1/d  below  0.5,  as  shown  by  Tables  II  and  III,  page  19.  However,  £q.  (7)  like- 
rise  requires  the  determination  of  that  value  of  n  which  gives  the  minimum  collaps¬ 
ing  pressure  and  hence  is  somewhat  cumbersome  and  indirect. 

Equation  (a)  presents  a  striking  contrast  to  equations  (6)  and  (7).  It  is 
very  simple  and  gives  the  collapsing  pressure  directly  without  the  use  of  n.  More— 
over,  Eq.  (a)  checks  Eq.  (6)  even  more  closely  than  does  Eq.  C*7),  and  can  replace 
(?-)  in  all  practical  calculations. 

Best  Available  Cop\ 
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